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ABSTRACT 

In 2013 April a new magnetar, SGR 1745—2900, was discovered as it entered an out¬ 
burst, at only 2.4 arcsec angular distance from the supermassive black hole at the 
Centre of the Milky Way, Sagittarius A*. SGR 1745—2900 has a surface dipolar mag¬ 
netic field of ~ 2 x 10 14 G, and it is the neutron star closest to a black hole ever 
observed. The new source was detected both in the radio and X-ray bands, with a 
peak X-ray luminosity Lx ~ 5 x 10 35 erg s _1 . Here we report on the long-term Chan¬ 
dra (25 observations) and XMM-Newton (8 observations) X-ray monitoring campaign 
of SGR 1745—2900 from the onset of the outburst in April 2013 until September 2014. 
This unprecedented dataset allows us to refine the timing properties of the source, as 
well as to study the outburst spectral evolution as a function of time and rotational 
phase. Our timing analysis confirms the increase in the spin period derivative by a fac¬ 
tor of ~2 around June 2013, and reveals that a further increase occurred between 2013 
Oct 30 and 2014 Feb 21. We find that the period derivative changed from 6.6xl0 -12 
s s” 1 to 3.3xl0~ n s s _1 in 1.5 yr. On the other hand, this magnetar shows a slow flux 
decay compared to other magnetars and a rather inefficient surface cooling. In par¬ 
ticular, starquake-induced crustal cooling models alone have difficulty in explaining 
the high luminosity of the source for the first ~200days of its outburst, and addi¬ 
tional heating of the star surface from currents flowing in a twisted magnetic bundle 
is probably playing an important role in the outburst evolution. 
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1 INTRODUCTION 

Among the large variety of Galactic neutron stars, 
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Table 1. Log of Chandra/ACIS-S and XMM-Newton/EPIC observations. Exposure times for the XMM-Newton observations are 
reported for the pn, MOS1 and M0S2 detectors respectively, and source net counts refer to the pn detector. 


Obs. ID 

MJD 

Start time (TT) 
(yyyy/mm/dd hh:mm:ss) 

End time (TT) 
(yyyy/mm/dd hh:mm:ss) 

Exposure time Source net counts 
(ks) (xlO 3 ) 

14702* 

56424.55 

2013/05/12 10:38:50 

2013/05/12 15:35:56 

13.7 

7.4 

15040** 

56437.63 

2013/05/25 11:38:37 

2013/05/25 18:50:50 

23.8 

3.5 

14703* 

56447.48 

2013/06/04 08:45:16 

2013/06/04 14:29:15 

16.8 

7.6 

15651** 

56448.99 

2013/06/05 21:32:38 

2013/06/06 01:50:11 

13.8 

1.9 

15654** 

56452.25 

2013/06/09 04:26:16 

2013/06/09 07:38:28 

9.0 

1.2 

14946* 

56475.41 

2013/07/02 06:57:56 

2013/07/02 12:46:18 

18.2 

7.1 

15041 

56500.36 

2013/07/27 01:27:17 

2013/07/27 15:53:25 

45.4 

15.7 

15042 

56516.25 

2013/08/11 22:57:58 

2013/08/12 13:07:47 

45.7 

14.4 

0724210201+ 

56535.19 

2013/08/30 20:30:39 

2013/08/31 12:28:26 

55.6/57.2/57.2 

39.7 

14945 

56535.55 

2013/08/31 10:12:46 

2013/08/31 16:28:32 

18.2 

5.3 

0700980101+ 

56545.37 

2013/09/10 03:18:13 

2013/09/10 14:15:07 

35.7/37.3/37.3 

24.9 

15043 

56549.30 

2013/09/14 00:04:52 

2013/09/14 14:19:20 

45.4 

12.5 

14944 

56555.42 

2013/09/20 07:02:56 

2013/09/20 13:18:10 

18.2 

5.0 

0724210501+ 

56558.15 

2013/09/22 21:33:13 

2013/09/23 09:26:52 

41.0/42.6/42.5 

26.5 

15044 

56570.01 

2013/10/04 17:24:48 

2013/10/05 07:01:03 

42.7 

10.9 

14943 

56582.78 

2013/10/17 15:41:05 

2013/10/17 21:43:58 

18.2 

4.5 

14704 

56588.62 

2013/10/23 08:54:30 

2013/10/23 20:43:44 

36.3 

8.7 

15045 

56593.91 

2013/10/28 14:31:14 

2013/10/29 05:01:24 

45.4 

10.6 

16508 

56709.77 

2014/02/21 11:37:48 

2014/02/22 01:25:55 

43.4 

6.8 

16211 

56730.71 

2014/03/14 10:18:27 

2014/03/14 23:45:34 

41.8 

6.2 

0690441801+ 

56750.72 

2014/04/03 05:23:24 

2014/04/04 05:07:01 

83.5/85.2/85.1 

34.3 

16212 

56751.40 

2014/04/04 02:26:27 

2014/04/04 16:49:26 

45.4 

6.2 

16213 

56775.41 

2014/04/28 02:45:05 

2014/04/28 17:13:57 

45.0 

5.8 

16214 

56797.31 

2014/05/20 00:19:11 

2014/05/20 14:49:18 

45.4 

5.4 

16210 

56811.24 

2014/06/03 02:59:23 

2014/06/03 08:40:34 

17.0 

1.9 

16597 

56842.98 

2014/07/04 20:48:12 

2014/07/05 02:21:32 

16.5 

1.6 

16215 

56855.22 

2014/07/16 22:43:52 

2014/07/17 11:49:38 

41.5 

3.8 

16216 

56871.43 

2014/08/02 03:31:41 

2014/08/02 17:09:53 

42.7 

3.6 

16217 

56899.43 

2014/08/30 04:50:12 

2014/08/30 15:45:44 

34.5 

2.8 

0743630201+ 

56900.02 

2014/08/30 19:37:28 

2014/08/31 05:02:43 

32.0/33.6/33.6 

9.2 

0743630301+ 

56901.02 

2014/08/31 20:40:57 

2014/09/01 04:09:34 

25.0/26.6/26.6 

7.8 

0743630401+ 

56927.94 

2014/09/27 17:47:50 

2014/09/28 03:05:37 

25.7/32.8/32.8 

7.7 

0743630501+ 

56929.12 

2014/09/28 21:19:11 

2014/09/29 08:21:11 

37.8/39.4/39.4 

11.7 


Observations already analysed by Rea et al. (2013a). An additional Chandra /HRC observation was carried out on 2013 
April 29. 

Chandra grating observations. 

+ XMM-Newton observations. 


magnetars constitute the most unpredictable class 
(Mereghetti 2008; Rea & Esposito 2011). They are 
isolated X-ray pulsars rotating at relatively long periods 
(P ~ 2— 12 s, with spin period derivatives P ~ 10 _15 -10~ 10 
s s -1 ), and their emission cannot be explained within the 
commonly accepted scenarios for rotation-powered pulsars. 
In fact, their X-ray luminosity (typically Lx ~ 10 33 -10 35 
erg s~ 4 ) generally exceeds the rotational energy loss rate 
and their temperatures are often higher than non-magnetic 
cooling models predict. It is now generally recognized 
that these sources are powered by the decay and the 
instability of their exceptionally high magnetic field (up to 
B ~ 10 14 — 10 15 G at the star surface), hence the name 
“magnetars” (Duncan & Thompson 1992; Thompson & 
Duncan 1993; Thompson, Lyutikov & Kulkarni 2002). 
Alternative scenarios such as accretion from a fossil disk 
surrounding the neutron star (Chatterjee, Hernquist & 
Narayan 2000; Alpar 2001) or quark-nova models (Ouyed, 
Leahy & Niebergal 2007a,b) have not been ruled out (see 
Turolla & Esposito 2013 for an overview). 


The persistent soft X-ray spectrum usually comprises 
both a thermal (blackbody, kT ~ 0.3-0.6 keV) and a non- 
thermal (power law, T ~ 2-4) components. The former is 
thought to originate from the star surface, whereas the latter 
likely comes from the reprocessing of thermal photons in a 
twisted magnetosphere through resonant cyclotron scatter¬ 
ing (Thompson et al. 2002; Nobili, Turolla & Zane 2008a,b; 
Rea et al. 2008; Zane et al. 2009). 

In addition to their persistent X-ray emission, magne¬ 
tars exhibit very peculiar bursts and flares (with luminosities 
reaching up to 10 46 erg s _1 and lasting from milliseconds to 
several minutes), as well as large enhancements of the persis¬ 
tent flux (outbursts), which can last years. These events may 
be accompanied or triggered by deformations/fractures of 
the neutron star crust (“stellar quakes”) and/or local/global 
rearrangements of the star magnetic field. 

In the past decade, extensive study of magnetars in out¬ 
burst has led to a number of unexpected discoveries which 
have changed our understanding of these objects. The de¬ 
tection of typical magnetar-like bursts and a powerful en- 
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Table 2. Timing solutions. Errors were evaluated at the lcr confidence level, scaling the uncertainties by the value of the rms (yx5) °f 
the respective fit to account for the presence of unfitted residuals. 


Solution 

Rea et al. (2013a) 

Kaspi et al. (2014) 

This work (Solution A) 

This work (Solution B) 

Epoch T 0 (MJD) 

56424.5509871 

56513.0 

56513.0 

56710.0 

Validity range (MJD) 

56411.6 - 56475.3 

56457 - 56519 

56500.1 - 56594.1 

56709.5 - 56929 

H(T 0 ) (s) 

3.7635537(2) 

3.76363824(13) 

3.76363799(7) 

3.7639772(12) 

P(To) 

6.61(4) x 10 -12 

1.385(15) x 10~ n 

1.360(6) x 10" 11 

3.27(7) x 10" 11 

P (s-i) 

4(3) x 10- 19 

3.9(6) x lO" 19 

3.7(2) x 10” 19 

(-1.8 ±0.8) x lO" 19 

u(T 0 ) (Hz) 

0.265706368(14) 

0.265700350(9) 

0.26570037(5) 

0.26567642(9) 

u(T 0 ) (Hz s- 1 ) 

-4.67(3) x lO" 13 

-9.77(10) x lO" 13 

-9.60(4) x 10 -13 

-2.31(5) x 10 -12 

v (Hz s" 2 ) 

-3(2) x 10 _2 ° 

-2.7(4) x 10" 2 ° 

-2.6(1) x 10“ 2 ° 

(1.3 ±0.6) x lO" 20 

rms residual 

0.15 s 

51 ms 

0.396 s 

1.0 ^iHz 

xl (d.o.f.) 

0.85 (5) 

1.27 (41) 

6.14 (44) 

0.66 (10) 



Time (MJD) 


f i t i£ 

f 

! i ,i 



56460 56500 56540 56580 56620 

Time (MJD) 
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Figure 1. Upper panel, temporal evolution of the spin frequency 
of SGR 1745—2900. The solution given by Rea et al. (2013a) is 
plotted as a green solid line. The blue and magenta solid lines 
show solutions A and B of this work, respectively. The blue dashed 
lines are the extrapolation of solution A over the time-span of so¬ 
lution B. The black line represents the fit over the whole time in¬ 
terval covered by observations (see text), while the vertical dashed 
lines refer to the times of the SGR-like short bursts detected by 
Swift /BAT (on 2013 April 25, June 7, and August 5). Central 
panel phase residuals with respect to solution A (labelled as 
<f>(t)A), evaluated over the time validity interval MJD 56500.1 
- 56594.1. Lower panel phase residuals with respect to solution 
B (labelled as vg(t)), evaluated over the time validity interval 
MJD 56709.5 - 56929. 


hancement of the persistent emission unveiled the existence 
of three low magnetic field (B < 4 x 10 13 G) magnetars (Rea 
et al. 2010, 2012, 2014; Scholz et al. 2012; Zhou et al. 2014). 
Recently, an absorption line at a phase-variable energy was 
discovered in the X-ray spectrum of the \ow-B magnetar 
SGR 0418+5729; this, if interpreted in terms of a proton cy¬ 
clotron feature, provides a direct estimate of the magnetic 
field strength close to the neutron star surface (Tiengo et al. 
2013). Finally, a sudden spin-down event, i.e. an anti-glitch, 
was observed for the first time in a magnetar (Archibald et 
al. 2013). 

The discovery of the magnetar SGR 1745—2900 dates 
back to 2013 April 24, when the Burst Alert Telescope 
(BAT) on board the Swift satellite detected a short hard 
X-ray burst at a position consistent with that of the super- 
massive black hole at the centre of the Milky Way, Sagit¬ 
tarius A* (hereafter Sgr A*). Follow-up observations with 
the Swift X-ray Telescope (XRT) enabled characterization 
of the 0.3-10 keV spectrum as an absorbed blackbody (with 
kT ~ 1 keV), and estimate a luminosity of ~ 3.9 x 10 35 
erg s^ 1 (for an assumed distance of 8.3 kpc; Kennea et al. 
2013a). The following day, a 94.5 ks observation performed 
with the Nuclear Spectroscopic Telescope Array ( NuSTAR) 
revealed 3.76 s pulsations from the XRT source (Mori et al. 
2013). This measurement was subsequently confirmed by a 
9.8 ks pointing on April 29 with the High Resolution Cam¬ 
era (HRC) onboard the Chandra satellite, which was able to 
single out the magnetar counterpart at only 2.4+0.3 arcsec 
from Sgr A*, confirming that the new source was actually 
responsible for the X-ray brightening observed in the Sgr A* 
region (Rea et al. 2013a). Follow-up observations in the 1.4- 
20 GHz band revealed the radio counterpart of the source 
and detected pulsations at the X-ray period (e.g. Eatough et 
al. 2013a; Shannon & Johnston 2013). The SGR-like bursts, 
the X-ray spectrum, and the surface dipolar magnetic field 
inferred from the measured spin period and spin-down rate, 
B p ~ 2 x 10 14 G, led to classify this source as a magnetar 
(Mori et al. 2013; Kennea et al. 2013; Rea et al. 2013a). 

SGR 1745—2900 holds the record as the closest neutron 
star to a supermassive black hole detected to date. The dis¬ 
persion measure DM = 1778 + 3 cm~ 3 pc is also the highest 
ever measured for a radio pulsar and is consistent with a 
source located within 10 pc of the Galactic Centre. Fur¬ 
thermore, its neutral hydrogen column density Ah ~ 10 23 
cm' 2 is characteristic of a location at the Galactic Centre 
(Baganoff et al. 2003). The angular separation of 2.4 ± 0.3 
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Figure 2. Pulse profiles of SGR 1745—2900 obtained from Chandra observations in the 0.3-10 keV energy range. Epoch increases from 
left to right, top to bottom. Two cycles are shown for clarity. 


arcsec from Sgr A* corresponds to a minimum physical sep¬ 
aration of 0.09 ± 0.02 pc (at a 95 per cent confidence level; 
Rea et al. 2013a) for an assumed distance of 8.3 kpc (see e.g. 
Genzel et al. 2010). Recent observations of the radio counter¬ 
part with the Very Long Baseline Array ( VLBA) succeeded 
in measuring its transverse velocity of 236±11 km s _1 at po¬ 
sition angle of 22 ± 2 deg East of North (Bower et al. 2015). 
If born within 1 pc of Sgr A*, the magnetar has a ^90 per 
cent probability of being in a bound orbit around the black 
hole, according to the numerical simulations of Rea et al. 
(2013a). 

SGR 1745—2900 has been monitored intensively in the 
X-ray and radio bands since its discovery. Three high-energy 
bursts were detected from a position consistent with that of 
the magnetar on 2013 June 7, August 5 by Swift/ BAT, and 
on September 20 by the INTErnational Gamma-Ray Astro¬ 
physics Laboratory (INTEGRAL ) (Barthelmy et al. 2013a,b; 
Kennea et al. 2013b,c; Mereghetti et al. 2013). Kaspi et al. 
(2014) reported timing and spectral analysis of NuSTAR 
and Swift/XRT data for the first ~4 months of the magne¬ 
tar activity (2013 April -August). Interestingly, an increase 
in the source spin-down rate by a factor ~ 2.6 was observed, 
possibly corresponding to the 2013 June burst. The source 
has been observed daily with Swift/XRT until 2014 Octo¬ 
ber, and its 2-10 keV flux has decayed steadily during this 
time interval (Lynch et al. 2015). 


Radio observations made possible a value of the rota¬ 
tional measure, RM = 66960 ± 50 rad m -2 , which implies 
a lower limit of ~ 8 mG for the strength of the magnetic 
field in the vicinity of Sgr A* (Eatough et al. 2013b). Ob¬ 
servations with the Green Bank Telescope showed that the 
source experienced a period of relatively stable 8.7-GHz flux 
between 2013 August and 2014 January and then entered a 
state characterized by a higher and more variable flux, until 
2014 July (Lynch et al. 2015). 

In this paper we report on the X-ray long-term moni¬ 
toring campaign of SGR 1745—2900 covering the first 1.5 yr 
of the outburst decay. In Section 2 we describe the Chandra 
and XMM-Newton observations and the data analysis. In 
Section 3 we discuss our results; conclusions follow in Sec¬ 
tion 4. 


2 OBSERVATIONS AND DATA ANALYSIS 

The Chandra X-ray Observatory observed SGR 1745—2900 
26 times between 2013 April 29 and 2014 August 30. The 
first observation was performed with the HRC to have the 
best spatial accuracy to localize the source in the crowded 
region of the Galactic Centre (Rea et al. 2013a). The re¬ 
maining observations were performed with the Advanced 
CCD Imaging Spectrometer (ACIS; Garmire et al. 2003) 
set in faint timed-exposure imaging mode with a 1/8 sub- 
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array (time resolution of 0.4 s), and in three cases with the 
High Energy Transmission Grating (HETG; Canizares et al. 
2005). The source was positioned on the back-illuminated 
S3 chip. Eight observations were carried out by the XMM- 
Newton satellite using the European Photon Imaging Cam¬ 
era (EPIC), with the pn (Striider et al. 2001) and the two 
MOS (Turner et al. 2001) CCD cameras operated in full- 
frame window mode (time resolution of 73.4 ms and 2.6 
s, respectively), with the medium optical blocking filter in 
front of them. A log of the X-ray observations is given in 
Table [Tj 

Chandra data were analysed following the standard 
analysis thread^] with the Chandra Interactive Analysis of 
Observations software package (CIAO, version 4.6; Fruscione 
et al. 2006). XMM-Newton data were processed using the 
Science Analysis Software (SAifi, version 13.5.0). For both 
Chandra and XMM-Newton data, we adopted the most re¬ 
cent calibration files available at the time the data reduction 
and analysis were performed. 

2.1 Timing analysis 

We extracted all Chandra and XMM-Newton/EPIC-pn 
source counts using a 1.5 and 15-arcsec circles, respectively, 
centred on the source position. Background counts were ex¬ 
tracted using a nearby circular region of the same size. We 
adopted the coordinates reported by Rea et al. (2013a), i.e. 
RA = 17 h 45 m 40M69, Dec = -29°00'29'/84 (J2000.0), to 
convert the photon arrival times to Solar System barycen- 
tre reference frame. The effects of the proper motion rela¬ 
tive to Sgr A* on the source position are negligible on the 
timescales considered for our analysis (best-fit parameters 
are 1.6 < /.i Q < 3.0 mas yr -1 and 5.7 < fj,s < 6.1 mas yr -1 
at a 95 per cent confidence level; Bower et al. 2015). 

To determine a timing solution valid over the time in¬ 
terval covered by the Chandra and XMM-Newton observa¬ 
tions (from 2013 April 29 to 2014 August 30; see Tabic 0, 
we first considered the timing solutions given by Rea et al. 
(2013a; using Chandra and Swift) and Kaspi et al. (2014; 
using NuSTAR and Swift). In the overlapping time inter¬ 
val, before 2013 June 14 (MJD 56457), both papers report 
a consistent timing solution (see first column in Tableland 
green solid line in the upper panel of Fig. [T]). Kaspi et al. 
(2014) then added more observations covering the interval 
between 2013 June 14 and August 15 (MJD 56457-56519), 
and observed a P roughly two times larger than the pre¬ 
vious value (see Table 0. The uncertainties on the Kaspi 
et al. (2014) solution formally ensure unambiguous phase 
connection until 2013 November 11 (MJD 56607), allowing 
us to extend this phase-coherent analysis with the data re¬ 
ported here, and follow the evolution of the pulse phases 
between 2013 July 27 and October 28 (MJD 56500-56594; 
after which we have a gap in our data coverage of about 115 
days; see Table 0. 

In this time interval, we measured the pulse phase at 
the fundamental frequency by dividing our observations in 
intervals of 10 ks and using the solution given by Kaspi 
et al. (2014) to determine univocally the number of cycles 

1 http://cxc. harvard.edu/ciao/threads/point like. 

2 http://xmm.esac.esa.int/sas/ 
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Figure 3. Temporal evolution of the pulsed fraction (see text 
for our definition). Uncertainties on the values were obtained 
by propagating the errors on the maximum and minimum count 
rates. Top panel : in the 0.3-10 keV band. Central panel: in the 

O. 3-3.5 band for the Chandra (black triangles) and XMM-Newton 
(red points) observations. Bottom panel: in the 3.5—10 band for 
the Chandra observations (black) and in the 3.5—5 (blue), 5—6.5 
(light blue) and 6.5—10 keV (green) ranges for the XMM-Newton 
observations. 

between the various observations. By fitting the measured 
pulse phases with a cubic function, we obtained the solution 
dubbed A in Tabled which shows only slight deviations with 
respect to the solution published by Kaspi et al. (2014), but 
extends until 2013 October 28 (MJD 56594). The period 
evolution implied by solution A is plotted with a blue solid 
line in the upper panel of Fig. \l\ Our Chandra and XMM- 
Newton observations allow us to confirm the change in the 

P, which increased by a factor of ^2 around 2013 June (i.e. 
about two months after the onset of the outburst in 2013 
April), and remained stable until at least 2013 October 28. 
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Figure 4. Pulse profiles of SGR1745—2900 obtained from XMM-Newton/FiPIC-pn observations. Two cycles are shown for clarity. Left: 
pulse profiles in the 0.3-10 keV energy band. Right: pulse profiles in the 0.3-3.5, 3.5-5, 5-6.5, 6.5-10 keV energy bands (from left to right, 
top to bottom) for the first four observations. Black, red, green and blue colors refer to the first, second, third and fourth observation, 
respectively. 


Formally, the accuracy of solution A should guarantee 
that phase coherence is not lost before 2014 March 3 (MJD 
56721), i.e. comprising the first observation available after 
the 115 day gap between MJD 56594.1 and MJD 56709.5. 
However, fitting the phases derived for that observation with 
solution A shows large residuals. These clearly indicate that 
solution A is not valid after the gap. To investigate this 
change in the spin evolution of the source, we measured the 
spin frequency for all the observations performed after the 
gap by fitting with a linear function the phases determined 
over time intervals with lengths ranging from 2 to 10 ks, 
depending on the source flux. The values for the frequen¬ 
cies we measured in this way after 2014 February 21 (MJD 
56709) are much smaller than those predicted by solution A 
(see blue dashed line in the upper panel of Fig. [l]) • 

To determine the spin evolution of the source after the 
115 day gap in the observations (i.e. from MJD 56709), we 
then fitted the values of the spin frequency with a quadratic 
function, obtaining the non-coherent solution B (see Ta¬ 
ble [2]), plotted in the upper panel of Fig. |T|with a magenta 
solid line. Unfortunately, this solution is not accurate enough 
to determine univocally the number of rotations between the 
various observations. Still, the trend followed by the spin fre¬ 
quency after the gap clearly deviates from that shown before 
2013 October 28 via solution A, indicating a further increase 
of the spin-down rate. In particular, the P has further in¬ 
creased by a factor of ~2.5, and the P is smaller than that 
measured by solution A, even if the large error prevents us 


from detecting a change in the sign of the P at high signifi¬ 
cance. 

The large changes in the timing properties of the source 
since the onset of the outburst are also shown by the fact 
that a quadratic function gives a poor fit for the spin fre¬ 
quency evolution over the whole time interval covered by the 
observations = 5.04 for 26 d.o.f.; see black solid line in 
the upper panel of Fig. [l]). 

Summarizing, we derive a phase coherent solution (so¬ 
lution A, see Tableland blue solid line in the upper panel 
of Fig. |T]) that is able to model the pulse phase evolution be¬ 
fore the 115 day observations gap starting at MJD 56600, 
and which is compatible with the solution given by Kaspi 
et al. (2014) for the partly overlapping interval MJD 56457 
- 56519. After the observation gap, solution A is no longer 
able to provide a good description of pulse phases, and we 
are only able to find a solution based on the analysis of the 
spin frequency evolution (solution B, see Tableland ma¬ 
genta solid line in the upper panel of Fig. [l]). 

We then use timing solution A (up to MJD 56594.1) 
and solution B (from MJD 56709.5 onwards) to fold all 
background-subtracted and exposure-corrected light curves 
at the neutron star spin period during the corresponding 
observation (see Figs [2] and [4]). This allows us to extract 
the temporal evolution of the pulsed fraction, defined as 
PF=[Max -Min]/[Max + Min] (Max and Min being the max¬ 
imum and the minimum count rate of the pulse profile, re¬ 
spectively). To investigate possible dependences on energy, 
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we calculate the pulsed fractions in the 0.3-3.5 and 3.5-10 
keV intervals for the Chandra observations and in the 0.3- 
3.5, 3.5-5, 5-6.5, 6.5-10 keV ranges for the XMM-Newton 
observations (see Fig. [3J. 

2.2 Spectral analysis of Chandra observations 

For all the Chandra observations, we extracted the source 
counts from a 1.5-arcsec radius circular region centred on 
SGR 1745—2900. This corresponds to an encircled energy 
fraction of ~ 85 per cent of the Chandra point spread func¬ 
tion (PSF) at 4.5 keV. A larger radius would have included 
too many counts from the Sgr A* PSF, overestimating the 
flux of SGR 1745—2900 with only a marginal increase of the 
encircled energy fraction (less than ~ 5 per cent). We ex¬ 
tracted the background counts using three different regions: 
an annulus (inner and outer radius of 14 and 20 arcsec, re¬ 
spectively), four 2-arcsec radius circles arranged in a square 
centred on the source, or a 1.5-arcsec radius circle centred 
on the source position in an archival Chandra/ ACIS-S ob¬ 
servation (i.e. when the magnetar was still in quiescence). 
For grating observations we considered instead a circle of 
radius 10 arcsec as far as possible from the grating arms but 
including part of the diffuse emission present in the Galactic 
Centre. 

For ‘non-grating’ observations, we created the source 
and background spectra, the associated redistribution ma¬ 
trix files and ancillary response files using the specextract 
tooS For the three grating observations, we analyzed only 
data obtained with the High Energy Grating (0.8-8 keV). 
In all cases SGR 1745—2900 was offset from the zeroth-order 
aim point, which was centered on the nominal Sgr A* coordi¬ 
nates (RA = 17 h 45 m 4ffi00, Dec = -29°00'28'. , 1 (J2000.0)). 
We extracted zeroth-order spectra with the tgextract tool 
and generated redistribution matrices and ancillary response 
files using mkgrmf and fullgarf, respectively. 

We grouped background-subtracted spectra to have at 
least 50 counts per energy bin, and fitted in the 0.3-8 keV 
energy band (0.8-8 keV for grating observations) with the 
XSPEC0 spectral fitting package (version 12.8.lg; Arnaud 
1996), using the y 2 statistics. The photoelectric absorption 
was described through the tbabs model with photoioniza¬ 
tion cross-sections from Verner et al. (1996) and chemical 
abundances from Wilms, Allen & McCray (2000). The small 
Chandra PSF ensures a negligible impact of the background 
at low energies and allows us to better constrain the value 
of the hydrogen column density towards the source. 

We estimated the impact of photon pile-up in the non¬ 
grating observations by fitting all the spectra individually. 
Given the pile-up fraction (up to ~ 30 per cent for the first 
observation as determined with WebPlMMS, version 4.7), we 
decided to correct for this effect using the pile-up model of 
Davis (2001), as implemented in xspec. According to ‘The 
Chandra ABC Guide to pile-up ’0 the only parameters al¬ 
lowed to vary were the grade-migration parameter (a?), and 
the fraction of events in the source extraction region within 

3 Ancillary response files are automatically corrected to account 
for continuous degradation in the ACIS CCD quantum efficiency. 

4 http://heasarc.gsfc.nasa.gov/xanadu/xspec/ 

5 http://cxc.harvard.edu/ciao/download/doc/pile-up_ abc.pdf 


the central, piled up, portion of the PSF. Including this com¬ 
ponent in the spectral modelling, the fits quality and the 
shape of the residuals improve substantially especially for 
the spectra of the first 12 observations (from obs ID 14702 
to 15045), when the flux is larger. We then compared our re¬ 
sults over the three different background extraction methods 
(see above) and found no significant differences in the pa¬ 
rameters, implying that our reported results do not depend 
significantly on the exact location of the selected background 
region. 

We fitted all non-grating spectra together, adopting 
four different models: a blackbody, a power law, the sum of 
two blackbodies, and a blackbody plus a power law. For all 
the models, we left all parameters free to vary. However, the 
hydrogen column density was found to be consistent with 
being constant within the errors^ among all observations 
and thus was tied to be the same. We then checked that the 
inclusion of the pile-up model in the joint fits did not alter 
the spectral parameters for the last 10 observations (from 
obs ID 16508 onwards), when the flux is lower, by fitting 
the corresponding spectra individually without the pile-up 
component. The values for the parameters are found to be 
consistent with being the same in all cases. 

A fit with an absorbed blackbody model yields y 2 = 
1.00 for 2282 degrees of freedom (d.o.f.), with a hydrogen 
column density Ah = 1.90(2) x 10 23 cm -2 , temperature in 
the 0.76-0.90 keV range, and emitting radius in the 1.2-2.5 
km interval. When an absorbed power law model is used 
(y 2 = 1.05 for 2282 d.o.f.), the photon index is within the 
range 4.2-4.9, much larger than what is usually observed 
for this class of sources (see Mereghetti 2008; Rea & Espos¬ 
ito 2011 for reviews). Moreover, a larger absorption value is 
obtained (Ah ~ 3 x 10 23 cm -2 ). The large values for the 
photon index and the absorption are likely not intrinsic to 
the source, but rather an artifact of the fitting process which 
tends to increase the absorption to compensate for the large 
flux at low energies defined by the power law. The addition 
of a second component to the blackbody, i.e. another black¬ 
body or a power law, is not statistically required (y 2 = 1.00 
for 2238 d.o.f. in both cases). We thus conclude that a sin¬ 
gle absorbed blackbody provides the best modelling of the 
source spectrum in the 0.3-8 keV energy range (see Table 
0 . 

Taking the absorbed blackbody as a baseline, we tried to 
model all the spectra tying either the radius or the temper¬ 
ature to be the same for all spectra. We found y 2 = 1.38 for 
2303 d.o.f. when the radii are tied, with Ah = 1.94(2) x 10 23 
cm -2 , Rbb = 1.991q o 5 km and temperatures in the 0.66- 
0.97 keV range. We found instead y 2 = 1.04 for 2303 d.o.f. 
when the temperatures are tied, with Ah = 1.89(2) x 10 23 
cm -2 , kTsB = 0.815(7) keV and radii spanning from ~1.1 
to ~3 km. The goodness of fit of the latter model improves 
considerably if the temperatures are left free to vary as well 
(E-test probability of ~ 2 x 10 -17 ; fitting the temperature 
evolution with a constant yields a poor y 2 = 2.8 for 24 d.o.f. 
in this case). We conclude that both the temperature and the 
size of the blackbody emitting region are varying. Zeroth- 
order spectral data of the three grating observations were 


6 Here, and in the following, uncertainties are quoted at the 90 
per cent confidence level, unless otherwise noted. 
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Table 3. Chandra spectral fitting results obtained with an absorbed blackbody model (y 2 = 1.00 for 2282 d.o.f.). The hydrogen column 
density was tied to be the same in all the observations, resulting in Ng = 1.90(2) X 10 23 cm -2 (photoionization cross-sections from 
Verner et al. (1996) and chemical abundances from Wilms et al. (2000)). The a is a parameter of the XSPEC pile-up model (see Davis 
2001 and “The Chandra ABC Guide to pile-up"). The pile-up model was not included when fitting the HETG/ACIS-S spectra (obs. ID: 
15040, 15651, 15654). The blackbody radius and luminosity are calculated assuming a source distance of 8.3 kpc (see e.g. Genzel et al. 
2010). Fluxes and luminosities were calculated after removing the pile-up model. All errors are quoted at a 90 per cent confidence level 
for a single parameter of interest (Ay 2 = 2.706). 


Obs. ID 

a 

kTgg 

Rbb 

1-10 keV absorbed flux 

1-10 keV luminosity 



(keV) 

(km) 

( 10 -12 erg cm -2 s _1 ) 

( 10 35 erg s” 1 ) 

14702 

0.47(6) 

0.87(2) 

2 c + 0.2 

z - D -o.i 

l6.5±S;g 

4.7 (3) 

15040 

- 

0.90(2) 

2.5(1) 

15.5±°;° 3 

4.7 (4) 

14703 

0 47+0-07 
_0.06 

0.84(2) 

2 . 6 ( 1 ) 

1 9 7+O.5 

1 Z -'- 0.6 

3.9 (3) 

15651 

- 

0.87(3) 

2.4(2) 

-,9 c+0.07 

3.8 (4) 

15654 

- 

0.88(4) 

2.4(2) 

12.4±g;g 5 

3.5 (4) 

14946 

0.43(8) 

0.82(2) 

2.5(1) 

10.518;* 

3.3(3) 

15041 

0 42+ 0 06 
u, ^ z -0.05 

0.83(1) 

9 99 + O.IO 
^•^_o.09 

Q q+0.2 
y *°—0.3 

2 q + 0.2 

z,y -0.4 

15042 

0.55(7) 

0.83(1) 

2.14(9) 

8.3(3) 

2 e + 0.2 

z -°-0.4 

14945 

0.4(1) 

0.85(2) 

1.9(1) 

7 «+0.3 
' -°-0.4 

2.3(2) 

15043 

0.51(8) 

0.82(1) 

2 09 +0 10 

7 9+0.2 

' ’ z —0.3 

2 4 +0.2 

z -^ -0.3 

14944 

0 . 6 ( 1 ) 

0.84(2) 

1.9(1) 

7.0(4) 

9 9+O.2 
z ‘ z —0.3 

15044 

O 

00 

1 + 
o o 

o o 

00 CO 

0.81(1) 

2 03 +0 ' 10 

^.uo_0 09 

6.4(2) 

9 1 +0.2 

z " 1 —0.3 

14943 

0.4(1) 

0.80(2) 

2 n + 0.2 
Z ' U - 0.1 

6 1 + 0.2 
d , 1 - 0.4 

2.0(3) 

14704 

0.5(1) 

0.80±g;g? 

2 . 0 ( 1 ) 

6 O” 1 ” 0 ' 2 
o,u -0.3 

2.0 ( 2 ) 

15045 

0.42(9) 

0.82l°;° 2 

1.88(9) 

5 9 ”*” 01 
o,y - 0.2 

1 9 + 0 - 1 
- 0.2 

16508 

0 . 6 ( 2 ) 

0.80l°;°? 

1 fi^ + 009 
i.oo _ 0 io 

0 Q + 0.1 

0.2 

1 Q + 01 
- 0.2 

16211 

0.3(2) 

0.79(2) 

1 64 +0 ' 10 
-L-o^-o.og 

3 6 +01 
o.O—o 2 

1.2 ( 2 ) 

16212 

0.4(2) 

0.80(2) 

1 5i + 010 
± - o± -0.09 

3.2(1) 

1 . 1 ( 1 ) 

16213 

0.3(2) 

0.79(2) 

1 49+ 008 
^•^-0.07 

3.1(1) 

1.1 ( 1 ) 

16214 

0.4(2) 

0.79(2) 

1 43+ 010 
-L-^O_0.09 

2 . 8 ( 1 ) 

1.0 ( 1 ) 

16210 

0.4(2) 

0.82(3) 

1.3(1) 

2.70(7) 

0.9(1) 

16597 

0.5(2) 

0.76(3) 

1 4 +°- 2 

2 20 + °' 04 
^.^u_o 05 

0.8 ( 1 ) 

16215 

0.4(2) 

0.80(2) 

-1 99+O.U9 

i.zz-o 08 

2.11(4) 

0.7 (1) 

16216 

0.3(2) 

0.76(2) 

1 34+ 010 
-L-^-o.og 

1.91(4) 

0.7(1) 

16217 

0.3(2) 

0.76(2) 

1.3(1) 

1.80(3) 

0.67(9) 


fitted together and indipendently with this model, without 
including the pile-up component and fixing Ng to that ob¬ 
tained in non-grating fit: 1.9 x 10 23 cm” 2 (see Tableland 

Fig. El). 

2.3 Spectral analysis of XMM-Newton 
observations 

For all the XMM-Newton observations, we extracted the 
source counts from a circular region of radius 15 arcsec cen¬ 
tred on the source PSF, and the background counts through 
the same circle at the same position in an archival (2011) 
XMM-Newton observation of the Galactic Centre (obs. ID 
0694640301), when the magnetar was not detected and no 
transient events were identified within the source PSF. We 
built the light curves for the source and background event 
files to visually inspect and filter for high particle back¬ 
ground flaring in the selected regions. We checked for the 
potential impact of pile-up using the epatplot task of SAS: 
the observed pattern distributions for both single and dou¬ 


ble events are consistent with the expected ones (at a 1 <t 
confidence level) for all the three cameras, proving that the 
XMM-Newton data are unaffected by pile-up. 

We restricted our spectral analysis to photons having 
FLAG = 0 and pattern ^ 4(12) for the pn (MOSs) data 
and created spectral redistribution matrices and ancillary 
response files. We co-added the spectral files of consecutive 
observations (obs. ID 0743630201-301 and 0743630401-501; 
see Table [TJ to improve the fit statistics and reduce the back¬ 
ground contamination. We then grouped the source spectral 
channels to have at least 200 counts per bin and fitted the 
spectra in the 2-12 keV range, given the high background 
contamination within the source PSF at lower energies. The 
spectral data extracted from the two MOS cameras gave 
values for the parameters and fluxes consistent with those 
obtained from the pn camera. To minimize the systematic 
errors introduced when using different instruments, we con¬ 
sidered only the pn data, which provide the spectra with the 
highest statistics. 

Due to the large PSF of XMM-Newton, it is not possible 
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Figure 5. Temporal evolution of the spectral parameters for the blackbody model and of the 1-10 keV absorbed flux of SGR 1745—2900 
from Chandra observations. 



Figure 6. Results of the phase-averaged spectral analysis for the XMM-Newton/EPIC-pn observations of SGR 1745—2900. Left panel. 
source spectra fitted together with an absorbed blackbody plus power law model in the 2-12 keV range and after removal of the Fe XXV 
and S XV lines (see text). E 2 X f(E) unfolded spectra together with the contributions of the two additive components and residuals 
(in units of standard deviations) are also shown. Right panel: source spectra fitted together with an absorbed 3-D resonant cyclotron 
scattering model in the 2—12 keV range and after removal of the Fe XXV and S XV lines (see text). Residuals (in units of standard 
deviations) are also shown. 


to completely remove the contamination of both the Galac¬ 
tic Centre soft X-ray diffuse emission and the emission lines 
from the supernova remnant Sgr A East, including in partic¬ 
ular the iron line (Fe XXV; rest energy of 6.7 keV) and the 
sulfur line (S XV; rest energy of 2.46 keV) (see e.g. Maeda et 
al. 2002; Sakano et al. 2004; Ponti et al. 2010, 2013; Heard 
& Warwick 2013). These features were clearly visible espe¬ 
cially in the spectra of the last observations, when the flux 


is lower, and prevented us from obtaining a good spectral 
modelling in xspec. We thus decided to discard the energy 
interval comprising the Fe XXV line (6.4-7.1 keV) for all the 
spectra, as well as that associated with the S XV line (2.3- 
2.7 keV) for the spectrum of the last observations (obs. ID 
0690441801, 0743630201-301, 0743630401-501), involving a 
loss of ~ 9 per cent in the total number of spectral bins. 

Based on the results of the Chandra spectral analysis, 
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Table 4. XMM-Newton/ EPIC-pn spectral fitting results obtained with an absorbed blackbody plus power law model (x 2 = 1.13 for 624 
d.o.f.) and an absorbed 3D resonant cyclotron scattering model (x 2 = 1-14 for 624 d.o.f.). flbulk denotes the bulk motion velocity of the 
charges in the magnetosphere and A<b is the twist angle. For both models the hydrogen column density was tied to be the same in all the 
observations, yielding N B = 1.86j)?'?? X 10 23 cm -2 for the former and N B = 1.86+??? x 10 23 cm -2 for the latter (photoionization cross- 
sections from Verner et al. (1996) and chemical abundances from Wilms et al. (2000)). The blackbody emitting radius and luminosity are 
calculated assuming a source distance of 8.3 kpc (see e.g. Genzel et al. 2010). Fluxes were determined with the CFLUX model in XSPEC. 
All errors are quoted at a 90 per cent confidence level for a single parameter of interest (Ax 2 = 2.706). 


BB+PL 







Obs. ID 

kT BB 

(keV) 

Rbb 

(km) 

r 

PL norm 
( 10 - 3 ) 

1 10 keV BB/PL abs flux 
( 10 -12 erg cm -2 s _1 ) 

1-10 keV BB/PL luminosity 
( 10 35 erg s _1 ) 

0724210201 

0.79(3) 

1.9(2) 

9 q + 0.5 

z.O-o.7 

4.511:? 

5.0(1) / 3.3(2) 

1 7 +O .2 / i n+0.4 
^'-O.S / 4 U —0.3 

0700980101 

0.78(3) 

2 . 1 ( 2 ) 

i 7 + 0.8 
—1.3 

< 6.8 

5.7(1) / 2.2(2) 

2.0(3) / 0.5(3) 

0724210501 

0.79(4) 

9 I +0.3 
z,± - 0.2 

9 q + 0.5 
^•- 0.6 

< 4.5 

5.8(1) / 1.8(2) 

2 Q + 0 - 1 / 0 3 + 0,3 
z,u — 0.2 / u,o - 0.2 

0690441801 

n 79+0.03 

u - tz -0.04 

1.6(3) 

2 6 +0 ' 5 
z>o - 0.8 

4 5 + 8 - 8 
4 '°-3.8 

1.9(1) / 2.1±0;| 

0.8(2) / 0.8(3) 

0743630201-301 

0.71(6) 

1.3(3) 

9 -1+0.7 
z, 1 -1.4 


1.2(1) / 1.7(2) 

0.5(2) / 0.4(3) 

0743630401-501 

0 07 + 0 - 10 
u -°‘-0.07 

1.2(5) 

2 0 +0 ' 4 
Z ’ U -0.7 

e-si 

0.7(1) / 2.0(4) 

0.3(2) / < 0.5 

NTZ 







Obs. ID 

kT 

(keV) 

Pbulk 

A 0 
(rad) 

NTZ norm 

(io- 1 ) 

1-10 keV abs flux 
( 10 -12 erg cm -2 s _1 ) 

1-10 keV luminosity 
( 10 35 erg s -1 ) 

0724210201 

0.85(2) 

n 79 +O.O 9 
u - ‘ z -0.40 

0 40"*” 0 04 
u *^ u -0.24 

i-62tS;E 

8.3(1) 

2.5(2) 

0700980101 

0-85t°;° 2 

0 70 + 0 - 04 
u -' u —0.34 

0 40"*"® 03 
U.^U-o.23 

i-58±g;ii 

7.9(1) 

2.3(2) 

0724210501 

0.84(2) 

0 . 6 ( 2 ) 

0 41+ 0 - 02 
u -^ 1 -0.25 

1.5(1) 

7.6(1) 

2.3(3) 

0690441801 

0 77+ 0 - 04 
u - 1 ‘ -0.06 

n k+0.3 
UiO -0.2 

0 42" 1 " 0 - 06 
u -^ z -0.25 

0 q4+ 010 
u - y4 —0.07 

4.0(1) 

1.3(2) 

0743630201-301 

0 76+ 0 - 07 
u - ' D -0.10 

> 0.2 

0.43l° 0 ;“ 

o.eil?:?? 

2.9(1) 

0.9(3) 

0743630401-501 

0 fi t s+ 0 07 
u.oo _ 0 24 

0 32” 1 ” 0 ' 11 

u.O^,_o 09 

0-60t°+ 7 

0 . 681? 22 

2.7(1) 

0.9(3) 


we fitted the data first with an absorbed blackbody model. 
The hydrogen column density was consistent with being con¬ 
stant at a 90 per cent confidence level among all observa¬ 
tions and was tied to be the same in the spectral fitting. 
We obtained xt = 2.2 for 636 d.o.f., with large residuals 
at high energies. The latter disappear if an absorbed power 
law component is added and the fit improves considerably 
(xt = 1-13 for 624 d.o.f.; see left panel of Fig. [UJ . A fit 
with a two-blackbody model is statistical acceptable as well 
(xt = 1.13 for 624 d.o.f.) and yields temperatures of ~ 2 — 4 
keV and emitting radii of ~ 0.04 — 0.12 km for the second 
blackbody. However, this model would be physically hard to 
justify, since it is unlikely that these large temperatures can 
be maintained on a neutron star surface for such a long time. 
As an alternative to these fits, we applied a 3D resonant cy¬ 
clotron scattering model (NTZ: Nobili et al. 2008a,b; Zane 
et al. 2009), obtaining xt = 1.14 for 624 d.o.f. (see right 
panel of Fig. ©. The hydrogen column densities and fluxes 
inferred both from the BB+PL and the NTZ models are 
consistent with each other within the errors (see Table HI). 
To test the robustness of our results, we compared the in¬ 
ferred parameters with those derived by fitting the spectra 
without filtering for the spectral channels and applying the 
VARABS model for the absorption, which allows the chemi¬ 
cal abundances of different elements to vary (only the sulfur 
and iron abundances were allowed to vary for the present 
purpose). We found consistent values over the two methods. 


We conclude that both models successfully reproduce 
the soft X-ray part of the SGR 1745—2900 spectra up to ~12 
keV, implying that, similar to other magnetars, the repro¬ 
cessing of the thermal emission by a dense, twisted magneto¬ 
sphere produces a non-thermal component. The power law 
detected by XMM-Newton is consistent with that observed 
by NuSTAR (Kaspi et al. 2014), and its very low contribu¬ 
tion below 8keV is consistent with its non-detection in our 
Chandra data. 

2.4 Pulse phase-resolved spectral (PPS) analysis 

To search for spectral variability as a function of rotational 
phase and time, we first extracted all the spectra of the 
Chandra observations selecting three pulse phase intervals 
(see Fig. [5]): peak (<(>=0.5-0.9), minimum (0=0.2-0.5), and 
secondary peak (0=0.9-1.2). We adopted the same extrac¬ 
tion regions and performed the same data analysis as for the 
phase-averaged spectroscopy. 

For each of the three different phase intervals, we fit¬ 
ted the spectra of all Chandra observations jointly in the 
0.3-8 keV energy band with an absorbed blackbody model 
and tying the hydrogen column density to be the same in 
all the observations (the pile-up model was included). Since 
the values of the column density are consistent with being 
the same at a 90 per cent confidence level (1.90(4) x 10 23 
cnW 2 , 1.82(4) x 10 23 cm -2 , and 1.83:??' ?® x 10 23 cm -2 for the 
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Figure 7. Evolution of the blackbody temperatures (left panel) 
and radii (right panel) for the peak (black points) and the min¬ 
imum (red points) of the pulse profile for the Chandra observa¬ 
tions. 


peak, the secondary peak and the minimum, respectively), 
we fixed /Vh to 1.9 x 10 23 cm -2 , i.e. to the best-fit value 
determined with the phase-averaged spectroscopy (see Ta¬ 
ble QJ[). We obtained a good fit in all cases, with xt — 1-04 
for 1005 d.o.f. for the peak, xt = 1.10 for 635 d.o.f. for the 
secondary peak, and xt = 0.99 for 713 d.o.f. for the pulse 
minimum. The fit residuals were not optimal for energies 
>6 — 7 keV for the peak spectra, due to the larger pile-up 
fraction. We extracted the source counts excluding the cen¬ 
tral piled up photons (within a radial distance of 0.7 arcsec 
from the source position), and repeated the analysis for the 
peak spectra: the residuals are now well shaped, and the 
inferred values for the spectral parameters did not change 
significantly. 

The temporal evolution of the blackbody temperature 
and radius for both the peak and the pulse minimum are 
shown in Fig. |T| No particular trend is observed for the in¬ 
ferred temperatures, whereas the size of the emitting region 
is systematically lower for the pulse minimum. This is con¬ 
sistent with a viewing geometry that allows us to observe 
the hot spot responsible for the thermal emission almost en¬ 
tirely at the peak of the pulse profile, and only for a small 
fraction at the minimum of the pulsation. 

The higher statistics of the XMM-Newton/EPIC-pn 
data allowed us to put more stringent constraints on the 
variations of the X-ray spectral parameters along the spin 
phase. We extracted the background-subtracted spectra in 


six different phase intervals for each observation, as shown 
in Fig. [8] We fitted all spectra with a BB+PL model, adopt¬ 
ing the same prescriptions used for the phase-averaged spec¬ 
troscopy in the filtering of the spectral channels. We tied the 
hydrogen column density and the power law photon indices 
to the best-fit values determined with the phase-averaged 
analysis (see Table 0. We obtained statistically acceptable 
results in all cases. The evolutions of the blackbody tem¬ 
perature and emitting radius as a function of the rotational 
phase for all the observations are shown in Fig. [8] Variability 
of both the parameters along the rotational phase is more 
significant during the first observation (a fit with a constant 
yields xt — 2.6 for 5 d.o.f. in both cases) than in the follow¬ 
ing observations (xt ^ 1.4 for 5 d.o.f. in all cases). 

To search for possible phase-dependent absorption fea¬ 
tures in the X-ray spectra of SGR 1745—2900 (similarly to 
the one detected in SGR 0418+5729; Tiengo et al. 2013), 
we produced images of energy versus phase for each of the 
eight EPIC-pn observations. We investigated different en¬ 
ergy and phase binnings. In Fig. [9] we show the image for 
the observation with the highest number of counts (obs. ID 
0724210201), produced by binning the source counts into 100 
phase bins and 100-eV wide energy channels. The spin pe¬ 
riod modulation is clearly visible, as well as the large photo¬ 
electric absorption below 2keV. For all observations we then 
divided these values first by the average number of counts in 
the same energy bin and then by the corresponding 0.3-10 
keV count rate in the same phase interval. No prominent 
features can be seen in any of the images. 


2.5 X-ray brightness radial profiles 

For all the Chandra observations, we used the Chandra Ray 
Tracer (chartQ; Carter et al. 2003) to simulate the best 
available PSF for SGR 1745—2900, setting the exposure time 
of each simulation equal to the exposure time of the cor¬ 
responding observation. For the input spectrum in CHART 
we employed the blackbody spectrum of Table 02 account¬ 
ing for the pile-up. We then projected the PSF rays on to 
the detector plane via the Model of AXAF Response to X- 
rays software (MARX0, version 4.5.0; Wise et al. 2003). We 
extracted the counts of both the simulated PSFs and the 
ACIS event files through 50 concentric annular regions cen¬ 
tred on the source position and extending from 1 to 30 pixels 
(1 ACIS-S pixel corresponds to 0.492 arcsec). We then gen¬ 
erated the X-ray brightness radial profiles and normalized 
the nominal one (plus a constant background) to match the 
observed one at a radial distance of 4 pixels, i.e. at a dis¬ 
tance at which pile-up effects are negligible. A plot of the 
observed and simulated surface brightness fluxes (in units of 
counts x pixel -2 ) versus radial distance from the position 
of SGR 1745—2900 is shown in Fig. |TD] for the observation 
with the highest number of counts (obs. ID 15041). 

Extended emission around SGR 1745—2900 is clearly 
detected in all the observations, and it is likely dominated 
by the intense Galactic Centre diffuse emission. A detailed 
analysis of the diffuse emission, including its spatial exten- 


7 http://cxc.cfa.harvard.edu/chart. 

8 http://space.mit.edu/CXC/MARX. 
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Figure 8. Evolution of the blackbody temperatures (left) and radii (right) as a function of the rotational phase for the XMM-Newton 
observations. Spectra of consecutive observations were coadded (obs. ID 0743630201-301 and 0743630401-501: see the two lower panels). 
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Figure 9. Energy versus phase image for the XMM-Newton observation with the highest number of counts (obs. ID 0724210201). The 
image was obtained by binning the EPIC-pn source counts into 100 phase bins and energy channels of 100 eV, to better visualize the 
shape of the pulse profile and its dependence on energy. 


sion and spectral properties, is beyond the scope of this pa¬ 
per, and will be published in a subsequent work. 


3 DISCUSSION 

3.1 Outburst evolution and comparison with 
other magnetars 

The past decade has seen a great success in detecting mag- 
netar outbursts, mainly thanks to the prompt response and 
monitoring of the Swift mission, and to the dedicated follow¬ 
up programs of Chandra, XMM-Newton, and more recently, 
NuSTAR. The detailed study of about ten outbursts has 
shown many common characteristics (see Rea & Esposito 
2011 for a review; see also Fig. ED, although the precise trig¬ 
gering mechanism of these outbursts, as well as the energy 
reservoir responsible for sustaining the emission over many 
months, remains uncertain. 

All the outbursts that have been monitored with suffi¬ 
cient detail are compatible with a rapid (<days) increase in 


luminosity up to a maximum of a few 10 35 erg s -1 and a 
thermally dominated X-ray spectrum which softens during 
the decay. In the case of SGR 0501+4516 and IE 1547—5408, 
a non-thermal component extending up to 100-200 keV ap¬ 
pears at the beginning of the outburst, and becomes unde¬ 
tectable after weeks/months (Rea et al. 2009; Bernardini et 
al. 2011; Kuiper et al. 2012). 

The initial behavior of the 2013 outburst decay of 
SGR 1745—2900 was compatible with those observed in 
other magnetars. The outburst peak, the thermal emission 
peaked at about 1 keV, the small radiating surface (about 
2km in radius), and the overall evolution in the first few 
months, were consistent with the behavior observed in other 
outbursts. However, after an additional year of X-ray mon¬ 
itoring, it became clear that the subsequent evolution of 
SGR 1745—2900 showed distinct characteristics. The source 
flux decay appears extremely slow: it is the first time that we 
observe a magnetar with a quiescent luminosity < 10 34 erg 
s' 1 remaining at a luminosity > 10 35 erg s _1 for more than 
one year, and with a temperature decreasing by less than 
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10% from the initial ~ 1 keV . A further interesting feature 
of this source is that the non-thermal component (as de¬ 
tected by XMM-Newton) persisted on a very long temporal 
baseline during the outburst evolution. The flux due to the 
power law component does not change significantly in time 
and, as a result, its fractional contribution to the total flux 
is larger at late times: ~ 520 d after the outburst onset, 
< 50 per cent of the 1 10 keV absorbed flux is due to the 
non-thermal component. 

We first modelled the decay empirically to gauge the 
characteristic decay timescales. We adopted three different 
functions to model the blackbody temperature, radius and 
1 -10 keV absorbed flux temporal evolutions (see Fig. [5]): (i) 
a linear model; (ii) a power law: /(f) = fo,PL x t -r ; (iii) 
an exponential: f(t) = fo,ex P x exp[—(t — to)/r], where to is 
the epoch of the first burst detected (which we fixed to 2013 
April 24 in all cases) and r is the e-folding time. 

The temporal evolution of the magnetar temperature 
is well represented by a linear model [xt = 0.7 for 23 
d.o.f.), with initial temperature IcTbb, o = 0.85(1) keV and 
slope (—1.77 ± 0.04) x 10 -4 . The hot spot shrinking is 
best modeled by an exponential (yf, = 0.8 for 23 d.o.f.). 
Best-fit parameters are r = 640 ± 62 d and initial radius 
Rbb,o = 2.60 ± 0.08 km. The shape of the flux decline ap¬ 
pears to change in time and in fact none of these models 
can accurately describe the magnetar flux overall decay. The 
flux decay during the first 100 days since the outburst on¬ 
set is well modelled by a linear plus exponential model with 
r = 37 ± 2 d (xt = 1.5 for 4 d.o.f.). After ~ 100 days, 
the best fitting model turns out to be an exponential with 
r = 253 ±5 d (xl = 1.4 for 15 d.o.f.). 

3.2 Crustal cooling modelling 

We applied the crustal cooling model (see e.g. Pons & Rea 
2012) to the data collected during the 1.5-yr outburst of 
SGR 1745—2900. Although this model was successful in ex¬ 
plaining several other magnetar outbursts (Rea et al. 2012, 
2013b), in this case we could not reproduce the very slow 
cooling and high luminosity observed for this source. We ran 
several models varying the total injected energy, the angu¬ 
lar size, and the depth of the region where the energy is 
released, but we could not find any set of parameters that 
fit the data. 

In the framework of the starquake model, the maxi¬ 
mum temperature reached in the region where the energy 
is released is limited by neutrino emission processes. This 
internal temperature determines the maximum surface tem¬ 
perature and therefore the luminosity at which the out¬ 
burst peaks during the first few days. For injected energies 
> 10 43 erg, there is no signihcant increase in the peak lumi¬ 
nosity because the crustal temperature saturates (at about 
3 — 5 x 10 9 K) due to the efficient neutrino processes. After 
reaching the maximum luminosity (between 1 hour and 1 
day depending on the depth and injection rate), the cooling 
curve tracks the thermal relaxation of the crust. Indepen¬ 
dent of the initial injected energy and surface temperature, 
the luminosity is expected to drop below 10 35 erg s -1 af¬ 
ter <20-30 days (see e.g. Fig. 1 in Pons & Rea 2012), due 
to neutrino emission processes in the crust (mainly plasmon 
decay, and probably neutrino synchrotron for magnetar held 
strengths). 



Figure 10. Radial profile of the surface brightness for both 
the ACIS-S image of SGR 1745—2900 (red dots) and the 
CHARt/marx PSF plus a constant background (black dots) for 
the observation with the highest number of counts (obs. ID 
15041). The simulated surface brightness has been normalized to 
match the observed one at 4 pixels (one ACIS-S pixel corresponds 
to 0.492 arcsec). Extended emission around SGR 1745—2900 is 
clearly detected in all the observations. 


In Fig. [12] (left panel, lower curves) we show an exam¬ 
ple of the expected cooling curve of a magnetar with the 
same characteristics of SGR 1745—2900 . We assume that a 
sudden large energy release, E ~ 10 45 erg, heats up a layer 
of the outer crust up to 3 x 10 9 K. We also assume that 
the event affects the entire magnetar surface, to create the 
most favorable scenario (the luminosity simply scales with 
the area of the emitting region), and that the layer where 
the energy is injected extends from an external boundary at 
Pout ~ 3 x 10 9 gem -3 , to an inner boundary at pin ~ 2 
and 4 x 10 10 gem -3 (we show these two cases in the two 
lower curves of Fig. 1121 left panel). It is clear that, even in 
this most favorable case, the high luminosities observed at 
late times are difficult to reconcile with any cooling model. In 
particular, injecting more energy or changing pout will only 
affect the peak luminosity during the first days or weeks. On 
the other hand, injecting energy deeper into the crust (i.e. 
at higher Pin), is expected to change the late time evolution 
only slightly. This can be seen by comparing the solid and 
dashed lines in the left panel of Fig. 1121 which correspond 
to pin = 2 and 4 xlO 10 gem -3 , respectively. 

For illustrative purposes, we also show the cooling 
curves obtained when plasmon and synchrotron neutrino 
processes are switched off (see the upper curves in the left 
panel of Fig. H- These provide a much closer match to the 
data; however, there is no clear reason why these neutrino 
processes should not operate in these conditions. This exam¬ 
ple is only meant to highlight the relevance of understanding 
neutrino processes in the crust, especially under the pres¬ 
ence of strong fields. Another possibility to fit the data is to 
tune the energy injection, which must be maintained dur¬ 
ing the first ~ 200 days, resulting in a higher luminosity at 
late times. If we assume that only a region 5 km in radius 
is affected (this is closer to the ~ 2 km emitting region ob¬ 
served), we need a continuous injection of at least ~ 10 44 
erg s -1 (per day) for about 200 days, which results in a to¬ 
tal energy of a few 10 46 erg. While this energy budget may 
not be unrealistic, a physical mechanism that can operate 



14 F. Coti Zelati et al. 



Time (days from burst activation) 


Figure 11. Flux decay of all magnetar outbursts monitored with 
imaging instruments. Fluxes are absorbed in the 1—10 keV energy 
range (adapted and updated from Rea &; Esposito 2011). 


for such a long timescale is not known. A possibility might 
be a continuous injection of energy to keep the surface at 
high temperatures for so long, although in this latter case 
we should possibly expect more SGR-like bursts during the 
first hundreds days. 


3.3 Bombardment by magnetospheric currents in 
a bundle 

In this section we discuss the possibility that the prolonged 
high luminosity of SGR1745—2900 is in part due to external 
particle bombardment as a consequence of the existence of a 
twisted magnetic field bundle. A valid alternative model to 
the crustal cooling scenario invokes the presence of magne¬ 
tospheric currents flowing along a gradually shrinking mag¬ 
netic bundle, and heating the surface from outside. Accord¬ 
ing to Beloborodov (2007, 2013), this bundle can untwist on 
different timescales: i) in the equatorial regions of the mag¬ 
netosphere, where the magnetic held reaches a few stellar 
radii, currents are dissipated after weeks or months, while 
ii) at higher latitudes (close to the poles), a bundle may 
untwist more slowly, possibly in one to ten years. Here, par¬ 
ticles can reach Lorentz factors of a few tens (Beloborodov 
2007). In this scenario, a quasi steady-state outflow of elec¬ 
trons and positrons is maintained thanks to magnetic pair 
production close to the surface. The non-negligible electric 
voltage along the magnetic held lines and the radiative force 
due to Compton scattering regulate the streams of positrons 
and electrons along the held line. 

The presence of a non-thermal component observed by 
NuSTAR (Mori et al. 2013; Kaspi et al. 2014), and confirmed 
also by our XMM-Newton observations on a much longer 
temporal baseline, is suggestive of a large density of mag¬ 
netospheric particles which boost thermal photons emitted 
from the surface via resonant Compton scattering, provid¬ 
ing the power law component. In this context, the observed 
~2 km size of the emitting blackbody is consistent with a 
relatively small j-bundle. In the scenario in which the out¬ 
burst evolution is dominated by an untwisting bundle and 
the poloidal magnetic held has a dipole geometry, the lumi¬ 


nosity is expected to decrease with the square of the black- 
body area ( At, = 4tyR% b ; Beloborodov 2007, 2009). A hatter 
dependence may arise from a more complex held geometry. 
In Fig. [12] we show the hts of the bolometric luminosity as 
a function of R% B with two different models, a quadratic 
function Lbol oc A% (black line; xt — 1-3 for 23 d.o.f.) and 
a power law Lbol oc Af (red line; = 0.8 for 23 d.o.f.). 
For the latter model we find a = 1.23(8). Interestingly, a 
similar relation was observed also for the outburst decay of 
SGR 0418+5729 (Rea et al. 2013b) and CXOU J1647-4552 
(An et al. 2013). 

In the following we will assess, using hrst order approx¬ 
imations, whether the particle density needed to keep the 
footprint of the bundle at a temperature of ~lkeV for the 
hrst hundreds of days after the outburst onset is consistent 
with the particle density in the bundle responsible for the 
non-thermal power law tail. The power of the infalling par¬ 
ticles is -EkinA) where -Ekin is the kinetic energy of a single 
particle at the surface and N is the total number of infalling 
particles per unit time. If this kinetic energy is transferred by 
the infalling particles to the footprint of the bundle, and pro¬ 
duces thermal luminosity from the footprint surface, then: 

Lx = AbcrT 4 = Efo n N = nVm e c 3 Ab, ( 1 ) 


where Ab is the area of the footprint surface, T is the spot 
temperature, n is the density of the infalling particles (as¬ 
sumed to be electrons and/or positrons, created by means of 
pair production), and F is the Lorentz factor. We calculated 
the density of the infalling particles by considering the ki¬ 
netic energy they need to heat the base of the bundle spot. 
For a given temperature, one can estimate n as 
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On the other hand, we can estimate the density of the 
particles responsible for the resonant Compton scattering 
which produces the X-ray tail as 
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where Jb = (c/47r)V x B is the conduction current, B is 
the local magnetic field, and r is the length-scale over which 
B varies (R„ ~ 10 6 cm is the star radius). In the magneto¬ 
sphere of a magnetar the real current is always very close 
to Jb and it is mostly conducted by e* pairs (Beloborodov 
2007). The abundance of pairs is accounted for by the multi¬ 
plicity factor M which is the ratio between the actual charge 
density (including pairs) and the minimum density needed 
to sustain J B ; the latter corresponds to a charge-separated 
flow in which the current is carried only by electrons (and 
ions). If the same charge population is responsible for both 
resonant Compton scattering and surface heating, the den¬ 
sities given by eqs. © and © should be equal. This implies 

B “(iy‘ Mr ’ 2 ^ 10 ‘(rSv) 4 ' < 4 > 


According to Beloborodov (2013), both the Lorentz fac¬ 
tor and the pair multiplicity change along the magnetic field 
lines, with typical values of M ~ 100 (i.e. efficient pair cre¬ 
ation), r ~ 10 in the largest magnetic field loops, and M ~ 1 
(i.e. charge-separated plasma), F ~ 1 in the inner part of the 
magnetosphere. The previous equality cannot be satisfied for 
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Figure 12. Left panel: crustal cooling curves attempting at modelling the luminosity decrease of SGR 1745—2900. Luminosities are 
bolometric and calculated assuming a distance of 8.3 kpc. For the neutrino-cooling on and off set of curves, the lower and upper curves 
are relative to pin = 2 and 4 X lO 10 gem -3 , respectively. Right panel: bolometric luminosity as a function of the square of the blackbody 
radius at infinity. Solid lines represent the fits with a quadratic function (black) and a generic power law (a = 1.23(8); red). 


a typical temperature of ~ 0.8 — 1 keV, unless the magnetic 
field changes over an exceedingly small length-scale, a few 
meters at most. It appears, therefore, very unlikely that a 
single flow can explain both surface heating and resonant 
up-scattering. 


4 CONCLUSIONS 

The spectacular angular resolution of Chandra and the large 
effective area of XMM-Newton, together with an intense 
monitoring of the Galactic Centre region, has allowed us 
to collect an unprecedented dataset covering the outburst of 
SGR 1745—2900, with very little background contamination 
(which can be very severe in this region of the Milky Way). 

The analysis of the evolution of the spin period allowed 
us to find three different timing solutions between 2013 April 
29 and 2014 August 30, which show that the source period 
derivative has changed at least twice, from 6.6xl0 -12 s s _1 
in 2013 April at the outburst onset, to 3.3xl0~ n s s _1 in 
2014 August. While the first P change could be related with 
the occurrence of an SGR-like burst (Kaspi et al. 2014), 
no burst has been detected from the source close in time 
to the second P variation (although we cannot exclude it 
was missed by current instruments). This further change in 
the rotational evolution of the source might be related with 
the timing anomaly observed in the radio band around the 
end of 2013 (Lynch et al. 2015), unfortunately during our 
observing gap. 

The 0.3-8 keV source spectrum is perfectly modelled by 
a single blackbody with temperature cooling from ~0.9 to 
0.75 keV in about 1.5 years. A faint non-thermal component 
is observed with XMM-Newton. It dominates the flux at 
energies >8keV at all the stages of the outburst decay, with 
a power law photon index ranging from ~ 1.7 to ~ 2.6. It 
is most probably due to resonant Compton scattering onto 
non-relativistic electrons in the magnetosphere. 

Modelling the outburst evolution with crustal cooling 


models has difficulty in explaining the high luminosity of this 
outburst and its extremely slow flux decay. If the outburst 
evolution is indeed due to crustal cooling, then magnetic 
energy injection needs to be continuous over at least the 
first ~200 days. 

The presence of a small twisted bundle sustaining cur¬ 
rents bombarding the surface region at the base of the bun¬ 
dle, and keeping the outburst luminosity so high, appears a 
viable scenario to explain this particular outburst. However, 
detailed numerical simulations are needed to confirm this 
possibility. 

This source is rather unique, given its proximity to 
Sgr A*. In particular, it has a > 90 per cent probability 
of being in a bound orbit around Sgr A* according to our 
previous N-body simulations (Rea et al. 2013a), and the re¬ 
cent estimates inferred from its proper motion (Bower et al. 
2015). We will continue monitoring the source with Chandra 
and XMM-Newton for the coming year. 
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